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C/3 ABSTRACT 

( ^-» ) Observations from the Heliospheric Imagers (HI- 1) on both the STEREO spacecraft have been 

analysed to search for bright low mass eclipsing binaries (EBs) and potential brown dwarf 
transits and to determine the radii of the companions. A total of 9 EB candidates have been 
found, ranging in brightness from V = 6.59 mag to V — 11.3 mag, where the radius of the 
companion appears to be less than QAR & , with a diverse range of host temperatures, from 
4074 K to 6925 K. Both components of one candidate, BD-07 3648, appear to be less than 
0AR & and this represents a particularly interesting system for further study. The shapes of the 
eclipses in some cases are not clear enough to be certain they are total and the corresponding 
radii found should therefore be considered as lower limits. The EBs reported in this paper 
have either been newly found by the present analysis, or previously reported to be eclipsing 
by our earlier STEREO/Hl-l results. One of the new objects has subsequently been confirmed 
using archival SuperWASP data. This study was made possible by using an improved matched 
filter extraction algorithm, which is described in this paper. 
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1 INTRODUCTION 

Observations of bright host stars of low mass eclipsing companions 
facilitates precise measurements of the mass and radius of the stars 
in these systems ( Jo hnson et al.||20TT) . There are, unfortunately, 
just 18 detached, double-lined eclipsing binary systems known that 
include a component with less than 0.8 M@ ( |Feiden & Cha boyer 
|2012[ >. The primary objective of this work was to search for pre- 
viously unknown low-mass EBs from STEREO/Hl-l data and to 
determine the radii of the companions. This could not be attempted 
in the earlier analysis of Wraight et al. ( 20 1 1 ) as it was only through 
the subsequent development of a matched filter algorithm (see Sec- 
tion |2.2[ ( that the difficulties of using the STEREO HI- IB data (see 
Section \2A\ could be overcome, which facilitates the extractions 
of low amplitude signals from the combined data, which together 
provide improved phase coverage across the variability cycle. Com- 
pared to our earlier analysis reported in Wraight et al. (201 1|>, up- 
dated flatfields had since become available (Bewsher et al.|2010) 
and an additional year of data from both satellites for each star 
had been taken, providing a baseline of about four and a half years 
compared to the two or three years in our previous study. Since 
the initial work of Wraight et al.|pQ"TT| was restricted to analysing 
each year of data separately using STEREO HI-1A data only, im- 
proved accuracy of eclipse depths and periods could be obtained. 
The analysis with the matched filter algorithm has allowed us to 
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identify additional candidates that had been missed by the previous 
search. In total, nine candidates are presented where the compan- 
ions appear to be less than 0.4 R@ and where the eclipses are most 
likely total and with a low risk of contamination from nearby stars, 
although the potential for the eclipses to be grazing remains in some 
cases. 

Studies of very low mass stars close to and below the boundary 
where the interior is fully convective are limited by a lack of bright 
eclipsing systems involving such stars (Feiden & Chaboyer 2012| >. 
Where these systems are bright enough, spectroscopic follow-up 
observations can unambiguously reveal the masses and radii of both 
components, allowing other parameters to be calculated with suffi- 
cient precision to challenge models of stellar formation and evolu- 
tion (|Shkomik et al.|200"8l|Cakirli, Ibaoglu & Sipahi|2012||Luhman| 
|2012| |Triaud et a l.|2012) . Indeed, there is an apparent tendency for 
current models to under-estimate the radii of low-mass stars, which 



suggests the presence of a mechanism inflating their radii ( Cakirli 



Ibaoglu & Dervisoglu|20T0l|Morales et al.|2010; Morales, Ribas & 



Jordi 201 1 1. It is therefore highly desirable to have a diverse sample 



of bright, low-mass detached EBs to compare against these models. 

The primary mission of NASA's twin STEREO spacecraft is 
to study the Sun's corona in three dimensions (Kaiser et al. [ 20081. 
Data from the Heliospheric Imager on the STEREO-Ahead space- 
craft (STEREO/Hl-lA; \Eyles et al.|2009} has previou sly been used 
to produce a catalogue of EBs (Wrai ght et al.|20TT) , to study the 
rotational periods of magnetic chemically peculiar stars (jWraight 
|et al.||20T2") and the pulsational periods of long period variables 
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Figure 1. STEREO/Hl-IA (left) and POSSI (right) views of HIP 247 and 
HIP 248, illustrating the difficulties of blending with the 70 arcsec pixel 
scale of STEREO/Hl-l, 



(Wraight et al.|2012b> , which was possible owing to the stability of 
the photometry and the accuracy of the calibration (Brown, Bew- 
sher & Eyles 2009, |Bewsher et al.|2010 Bewsher, Brown & Eyles 
2012). Several new EBs from |Wraig ht et al.|p0Tl) were suspected 
to be low mass and confirmation of this provided the motivation 
to search through data from both the Ahead STEREO spacecraft 
(STEREO HI-1A) and the Behind spacecraft (STEREO HI- IB) us- 
ing a custom matched filter algorithm to extract weak variable sig- 
nals. 

This paper describes the matched filter algorithm used to carry 
out the large-scale analysis of STEREO/HI- 1 data in Sectionp] The 
candidate low mass EBs are presented in Section[3] In Section|4]the 
sample is discussed in the context of current research in this field. 

This sample of low mass EB candidates represents an example 
of one limited study of STEREO data. There are numerous possi- 
bilities for this data and the prospects for the legacy of STEREO 
observations are huge, with a large coverage on the sky of bright 
stars, regardless of spectral type. 



2 DATA ANALYSIS 



Data from the Heliospheric Imagers (STEREO/Hl-l) (Eyle s et al.| 
|2009| > on both STEREO spacecraft was used for the analysis. A 
matched filter algorithm (see Section [2~2] l was developed to search 
through the lightcurves of all of the stars in the STEREO database. 
The reason a matched filter was used was to allow a good signal 
to be obtained from a variety of variability types (Jenkins, Doyle 



& Cullers 1996, Weldrake & Sackett 2005, Enoch 2009; Enoch et 



aLpn p ress 1. However, due to blending with nearby stars, the re- 



gion close to the Galactic Centre was excluded from the search. 
The pixel scale of the STEREO/HI- 1 imagers is 70 arcsec per pixel 
( |Eyles et al.|2009| > and close to the Galactic Centre, the high den- 
sity of stars on the sky results in the eclipse depth to be less than the 
real depth because of blending in the same and/or from neighbour- 
ing pixels. To illustrate this, Figure [T] shows HIP 247 and HIP 248, 
as seen by STEREO/Hl-l A and compared with an image of the 
same stars from the POSSI survey. It is not possible in this par- 
ticular case to determine which of the two stars is the source of 
the observed shallow eclipses ( Wraig hTet al.|20lT} from STEREO 
data alone and several other candidates were similarly found to be 
unresolvable, necessitating follow-up measurements to confirm the 
companions' status as low mass EBs (most notably BD+03 263p / 
BD+03 263s and BD+03 2482 / BD+03 2483A). The matched filter 
algorithm we have used was developed to search for both general 
variability and detached and semi-detached EBs in stars brighter 
than R = 1 1th magnitude. 



Where available, archival SuperWASP ( TPollacco et al. |2006| > 
observations of these candidates were checked, to help determine 
whether any variability was due to a background EB, as additional 
confirmation of the primary eclipse depth and to clarify whether 
or not secondary eclipses were present. Whilst the pixel scale of 
the Heliospheric Imagers may result in blending in several of the 
candidates, the eclipse depths and shapes still suggest the presence 
of a low mass companion. Additional observations will be required 
to unambiguously rule out the possibility of grazing eclipses. Fur- 
thermore, radial velocity measurements are required to confirm the 
low-mass nature of the companions. We do however note that a sig- 
nificant amount of radial velocity data already exists for HD 75767 
l |Griffin|2006) and HD 213597 jNordstroem et al.|1997| >, both sets 
of data supporting our suggestion that these are low mass binary 
systems. 

2.1 Characteristics of STEREO data 



The STEREO Heliospheric Imagers (Eyles et al. 2009) and their 
calibration (Brown et al. 2009 ; Bewshe r et al.|2010||2012) are sum- 
marised here. The STEREO/Hl-l imagers have a field of view of 20 
degrees by 20 degrees, which is imaged onto a 2048 x 2048 element 
CCD, then binned on-board 2x2 to produce the final 1024x 1024 
images with 70x70 arcsec pixels. The spectral response of the filter 
is peaked between 630 nm and 730 nm (closely resembling a con- 
ventional i?-band filter), although it does allow some light through 
close to 400 nm and 950 nm (Eyles et al. 2009). Stars of 3rd mag- 
nitude or brighter were excluded from the survey, whilst stars down 
to R = 12th magnitude from the NOMAD catalogue (Zach arias et| 
|al.|2004| > were used as the basis for selection. Over the course of an 
orbit, almost 900,000 stars of 12th magnitude and brighter are im- 
aged within 10 degrees of the ecliptic plane. Since the HI-1 imager 
field of view is centred only 14 degrees away from the centre of the 
Sun, data on the sunward side of the field of view occasionally can 
be affected by solar activity and this region was masked out from 
the analysis. 

The STEREO-Ahead spacecraft is in an Earth-leading orbit 
with a semi-major axis of 0.95 AU, whilst the STEREO-Behind 
spacecraft is in an Earth trailing orbit with a semi-major axis of 
1.05 AU ( jKaiser et al.|2008} . This results in stars remaining in the 
field of view of the STEREO/Hl-l A imager for just over 19 days 
and in the field of view of the STEREO/Hl-IB imager for just over 
22 days, with new observations being recorded every 40 minutes. 
With this combination of regular photometry repeated every orbit, 
the STEREO mission is therefore very well suited to a wide vari- 
ety of stellar variability studies requiring continuous observations 
over many weeks (Wraight et al. 201 11. For this analysis, four such 
passes of a region of the sky from each satellite were available for 
most stars, resulting in many thousands of data points spanning up 
to nearly four and a half years. With more than twice the amount 
of data and nearly double the time-span used for the analysis of 
|Wraight et alT| l |20TT] >, it was possible to extract a larger sample of 
low mass EB candidates. The tolerance of the matched filter al- 
gorithm to artefacts in the data also allowed the complete set of 
data for each star to be analysed together, whereas those used in 
Wraight et al. (2011) were restricted to separately using data from 
each year's observations in isolation to each other. Improved flat- 
fields also became available (Bewsher et al ,|2010} , improving the 
efficacy of the data and minimising the amount of detrending of the 
data that was necessary. 

The STEREO/Hl-IB data often exhibits sudden decreases 
in amplitude of all signals at a given time (Figure [2}. These 
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are attributed to errors in the background subtraction due to 
pointing changes resulting from micrometeorite impacts, as the 
STEREOfHl-lB imager is facing the direction of travel along the 
spacecraft's orbit, recently confirmed by |Davis et al.|(2012| >. The 
severity of events can range from a single observation point to days 
of erratic behaviour. Often there will be a marked decrease for some 
hours and then a return to normal. In spite of the potential for mis- 
classification of such signals as eclipses, they can normally be eas- 
ily identified by visual examination, as the shape and duration of 
eclipses in a real EB do not change, beyond the phase coverage 
and accuracy of the observations, whereas these artefacts are errat- 
ically shaped and, importantly, aperiodic. However, such effects are 
relatively uncommon in STEREOfHl-lA data, and also less severe 
when they are present, making these observations much more stable 
than those of STEREO/H1-IB. The matched filter algorithm devel- 
oped for the analysis is not readily confused by aperiodic changes, 
since the signal strength found does not usually change for differ- 
ent periods. However, occasional extreme events cause problems 
for the polynomial detrending used to clean the lightcurve before 
being passed to this algorithm (see Section |272| below). 



2.2 Methods used in the analysis 

Before being passed to a custom matched filter algorithm for analy- 
sis, all lightcurves were passed through a culling routine to remove 
some of the artefacts caused by Mercury and Venus passing through 
the field of view as well as some of the more extreme pointing- 
related events, although as it is intended to leave some variabil- 
ity untouched, such as from eclipses, this is not wholly effective. 
Points more than 4a away from the weighted mean magnitude, for 
each epoch, were excluded in this process. Polynomial detrending 
was then carried out by a 4th order polynomial in order to remove 
residual trends from the flat fielding or other artefacts. This also 
removes some variability that is long compared to the length of an 
epoch, thus periods longer than about half the length of an epoch 
are not expected to be reliable, however long period eclipses are not 
removed in this process. 

The matched filter algorithm analyses a lightcurve in several 
stages, outlined below, and works by building model lightcurves 
and finding the minimum least-squared error of a model compared 
to the actual lightcurve. The model lightcurves consist of data 
points with the same time and errors as the real lightcurve. The 
sequence of the processing steps is as follows: 

(i) Determine a best-fitting period. A sinusoidal shape is used at 
this stage. A periodogram is produced during this process. 

(ii) Fine-tune the period. A precise period is the most important 
step in the process and without this the other characteristics will 
not be reliably determined. 

(iii) Determine the amplitude of variability. This needs to be 
done before and after determining the best-fitting shape. 

(iv) Determine the shape of the variability. A selection of shapes 
consisting of sinusoidal variability with different harmonic signals 
overlaid are used in addition to shapes based upon box-like total 
eclipses, v-shaped eclipses and a composite of total eclipses with 
wide ingress and egress phases. The amplitude is recalculated after 
this stage and, if an eclipsing model was best-fitting, the duration 
and depth of eclipses. Figure[3]shows the different eclipsing models 
as applied to EBs either from the present sample or from previously 
known sources (it does not produce a v-shaped model for any EB 
in the sample). 



(v) If the best-fitting amplitude is zero this means that the model 
finds the star to be constant. The period will still be returned. 

(vi) If an eclipsing model was best-fitting and the amplitude 
non-zero, eccentricity and amplitude of secondary eclipses is 
checked for. Higher harmonics of the period are also checked. 

This matched filter algorithm is processor-intensive and, de- 
pending on the number of data points, may take from 5 to 10 min- 
utes per star to process for a period range from 0.1 to 3.5 days. 
This has a resolution of 0.005 days initially and then is fine-tuned 
to 0.00005 days in a narrow search around the best-fitting period. 
Higher frequencies were excluded to avoid the Nyquist frequency 
of about 15.625 days -1 (0.064 days). This was the period range 
used in the search, however the program showed during testing a 
clear preference for finding the half-period harmonic of EBs, to the 
extent that it was specifically programmed to check at double and 
triple the periods initially found to see if these produced a better fit. 
It thus potentially checked periods up to 10.5 days for eclipses. 

A table of statistics was produced after analysing each field 
of view. These tables could then be statistically analysed to 
try to identify low mass EB candidates. Visual examination of 
the periodograms produced by the matched filter, along with 
lightcurves phase-folded on the best-fitting period and the de- 
trended lightcurves themselves, was then undertaken to produce a 
shortlist of signals for further examination. 

After the initial list of candidates had been compiled, the sur- 
rounding stars in the vicinity of low mass EB candidates were 
checked to assess whether there was likely to be any contamination 
from other nearby known bright EBs. This ruled out many candi- 
dates. Those remaining were cross-checked in the literature to see 
if they were already known prior to any analysis with STEREO. 
Archival SuperWASP ( Poll acco et al.|2 006 1 data was available for 
some candidates and this was further checked to confirm both the 
source and depth of eclipses. A more detailed study of the known 
parameters for each star was carried out using information available 
in the literature, which showed one candidate to be likely to be a 
giant and the companion therefore too large to be a low mass EB. A 
handful of the remaining candidates were then set aside for poten- 
tial follow-up but excluded from presentation here, as the source of 
eclipses needed to be more reliably determined or the eclipses ap- 
peared more likely to be grazing. These candidates were subjected 
to the same analysis as the nine that were left, however the analysis 
is of more limited use when the shape and depth of eclipses is not 
reliably determined (in particular, two apparently promising new 
candidate EBs, HD 18731 and HD 73470, appeared more likely to 
be grazing). 



3 LOW MASS EB CANDIDATES 
3.1 Lightcurves and data 

Stars are listed in Table [T] in order of increasing radius of the 
companion. Errors on the eclipse depths have been determined as 
the standard deviation of the out-of-eclipse lightcurves. The pe- 
riod analysis was done in PeranscQ The Modified Julian Date 
(MJD-2400000) of mid-eclipse is given for each candidate. Figure 
[4] shows the lightcurves of the candidates phase-folded on the pe- 
riods given and Figure [5] shows close-ups of the eclipses, to help 
clarify which are most likely to be total, since for grazing eclipses 

1 http://www.peranso.com 
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Figure 2. One block of data from STEREO/HI- IB of the bright low-mass EB HD 15161 . The eclipses are labelled but all the other sharp decreases in brightness 
are believed to be due to de-pointing events resulting from micrometeorite impacts. The smoother changes in brightness are a systematic effect caused by the 
polynomial detrending. The error bars shown are the 1 o errors on the individual data points. A block of data such as this would typically be completely 
excluded when doing an individual analysis, however as the eclipses are clearly evident a limited section can be considered for inclusion. 
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Figure 3. The different eclipse shapes utilised by the matched filter algorithm described in Section |272| The EBs are HD 213597 (box type, left), V818 Tau 
(v-shaped, middle, note that secondary eclipses are just visible) and HD 222891 (composite with one-tenth of the eclipse in totality, right). 



the radii found will be lower limits. Table [2] gives the times when 
the stars in Table[T]will next be viewed by STEREO/HI- 1 , to facili- 
tate observers who may wish to make observations simultaneously. 



3.2 Determination of radii 

One of the most crucial pieces of information which can be pre- 
cisely extracted from an EB system is the radius of the companion 
star, given the knowledge of the eclipse depth and host star radius. 
The radii of the stars in our sample were estimated from the stellar 
temperature and luminosity, using the Stefan-Boltzmann relation, 
shown in Equation[TJ where R is radius, L is luminosity and T e ff is 
effective temperature. It was assumed in all cases that the observed 
eclipses are total (thus the radii derived are lower limits) and the ra- 
dius of the companion then follows geometrically from the radius 
of the primary. 



A' 

r7 



'effO 
T eff* 



(1) 



Only one of the stars in our sample has previously been anal- 
ysed spectroscopically, HD 75767, therefore it was decided to de- 
termine r e ff photometrically. HD 75767 has been estimated to have 
a T e ff of 5823 K and a bolometric absolute magnitude of 4.60 mag 



in |Mishenina et al.| ( |2009| >. For all of our stars Johnson (B and 
V from SlMBAD), 2MASS (J, H and K s from |Skrutskie et al.| 
|2006| > and Tycho (B T and V T from |H0g et al.||2000|> photometry 
were available. The photometric calibrations by 'Casagrand e et al.| 
( |2010fr were then used to estimate T e ff from various colour indices. 
Where available Stromgren photometry from Hauck & Mermilliod 



od 



( 1998 1 was also used with calibrations by Casagrande et al. (2010 
Moon & Dworetsky ( 1985 1, Napiwotzki, Schoenberner & Wenske 
( |1993t , |Balona| ( |1994| >, |Ribas et al.| ( |1997| >, and |Castelli, Gratton 
|& Kurucz {1997). The effective temperature finally adopted for 



each star, listed in Table [T] is the median value obtained from each 
colour whilst the uncertainty given is the standard deviation from 
the mean. 

To derive a precise value of the stellar luminosity it is crucial 
to have a precise measurement of the stellar distance. We have re- 
stricted our luminosity estimates to only those available in the new 
reduction of the HIPPARCOS catalog, presented by jvan Leeuwen| 
l |2007) . For these stars the stellar luminosity was derived from the 
V Johnson magnitude using the bolometric correction by Balona 
( |1994) , the effective temperature previously obtained, and redden- 
ing by Amores & Lepine ( 2005 ). The uncertainty on the luminos- 
ity was derived taking into account the uncertainties on the distance 
given by |van Leeu wen ( 2007 1, an uncertainty in the bolometric cor- 
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rection of about 0.07 mag, and a reddening uncertainty of 0.01 mag 
l |Fossati et al.|2008| >. 

For those stars not present in the new reduction of the HIP- 
PARCOS catalog it was decided just to give a rough estimate of the 
stellar radius on the basis of the effective temperature derived from 
the photometry (stars marked with a * in Table[T] For this operation 
the values used are those given by Gray 1992 1. As the stellar radius 
was just a rough estimate, for these stars no uncertainty was de- 
rived, although an error bar of at least 50% is plausible. Therefore 
these values have to be taken with caution. 



Table 1. Comparison of low mass EB candidates observed by STEREO/Hl-l in order of increasing radius of the companion. The spectral type is as given in SlMBAD but where this has not been determined it is 
listed as NA. The eclipse depth of each is given in mmag with an error corresponding to the standard deviation of the out-of-transit lightcurve, as found from the phase-folded lightcurves shown in Figure [4] The 
Modified Julian Date (MJD-2400000) of the mid-eclipse time is given from the period analysis, which also provides the periods and their errors. The T c ff and its error for each star was determined photometrically 
as described in Section |3~2] with the median value given here and the error being the standard deviation. Radii marked with a * are estimated on the basis of T c ff rather than calculated from photometry as shown in 
S ection|3.2jNote that the companion radii are lower limits and assume that the eclipses are total. In the last column, each candidate is flagged as to whether it is previously known to be eclipsing from the analysis 
of |Wraight efal^TlTT] . 



Identity 


RA 


DEC 


Mag 


Sp type 


Depth 


MJD 


Period 




^host 


Krnmn 


Known 




(deg) 


(deg) 


(V) 




(mmag) 




(days) 


(K) 


(«©) 


(R s ) 




HD23765 


57.1205 


21.7975 


9.53 


F8 


49 ±17 


54213.264019 


1.6865 ±0.0004 


5746 ± 142 


0.98* 


^ 0.22 


no 


HD287039 


71.5819 


12.7448 


9.88 


F8 


75 ±32 


54224.548593 


2.2111 ±0.0004 


5310±145 


0.84* 


> 0.23 


no 


HD89849 


155.541 


6.21829 


9.15 


F8 


54 ±23 


54305.226101 


3.0781 ±0.0017 


5789 ±83 


1.0* 


> 0.23 


yes 


BD-07 3648 


203.727 


-8.44247 


11.3 


NA 


472± 131 


54353.170536 


2.5265 ±0.0004 


4074 ±918 


0.35 ±0.09 


> 0.24 ±0.36 


no 


HD75767 


133.068 


8.06293 


6.59 


G0+M3 


103 ±13 


54285.947481 


10.2478 ± 0.0043 


5649 ±245 


1.02 ±0.09 


> 0.33 ±0.12 


no 


HD 198044 


312.091 


-22.7407 


7.2 


F7V 


53±12 


54450.509489 


5.15195 ±0.0015 


61 16 ±236 


1.46±0.11 


5s 0.34 ±0.1 3 


yes 


HD205403 


323.766 


-3.7349 


8.02 


F5 


57 ±14 


54123.634818 


2.4449 ±0.0005 


6771 ±136 


1.46* 


> 0.35 


yes 


HD2 13597 


338.136 


1.58245 


7.81 


F0 


33 ±13 


54138.602903 


2.4238 ±0.0007 


6925 ±112 


1.96 ±0.06 


^ 0.36±0.12 


yes 


HD222891 


356.162 


-8.84879 


8.07 


F8 


39 ±12 


54150.650434 


1.59495 ±0.0004 


6409±112 


1.80 ±0.06 


^ 0.36±0.11 


yes 



Table 2. The dates shown (in MJD-2400000) in this table are when the stars will next be observed by STEREO/HI- 1 , until the two STEREO spacecraft have passed behind the Sun. The times take into account 
the mask that is routinely applied to exclude data likely to be affected by solar activity. As the Sun is going to be at maximum during these observations, it is not expected that observations will be reliable within 
the region of the CCDs excluded by the mask. The designation NA indicates that data will not be available for that star on that occasion, in particular HD 205403 is currently just outside the field of view of the 
STEREO/HI- IB imager and will not be visible during the coming orbits for that satellite. Observation periods that may be occurring when the relevant spacecraft are near to losing communication due to being 
behind the Sun are marked with a * and may not be reported, or suffer from a reduced cadence and missing data. 



Identity 



start and end dates of upcoming STEREO/HI- 1 observations (MJD-2400000) 



STEREO/UI-IA 



STEREO/HI-IB 



STEREO /HI-1 A 



STEREOim-lB. 



STEREO '/HI-l A 



STEREO '/HI- IB 



HD 23765 56280.01953-56293.18359 56202.00651 -56217.03385 56624.57552-56637.73958 56591.04297-56606.07031 

HD 287039 56291.18359-56310.37891 56239.84375-56259.09375 56635.68359-56654.87891 56642.99219-56662.24219 

HD 75767 56349.59766-56368.34766 56279.20508-56302.45508 56694.10677 -56712.85677 56668.17578-56691.42578 

HD 89849 56370.43750-56382.88281 56312.36328-56327.25000 56714.96615 -56727.41146 56701.26563 -56716.15234 

BD-07 3648 56417.92969-56431.03906 56368.05469-56382.02734 56762.46615 -56775.57552 56756.98438-56770.95703 

HD 198044 56172.34766-56185.67969 56473.01367-56486.29102 56516.87500-56530.20703 56861.95703 -56875.23438 

HD 205403 56187.16016-56204.74609 NA 56531.19922-56548.78516 NA 

HD 213597 56202.41016-56221.57813 56496.00391 -56514.67188 56546.53255 -56565.70052 56884.10156-56902.76953 

HD 222891 56214.27344-56225.60547 56124.37891 -56141.98828 56558.38672-56569.71875 56513.87891 -56531.48828 



56969.13151 -56982.29557 
56980.18359-56999.37891 
57038.61589 - 57057.36589 
57059.49479 - 57071.94010 
NA 

56861.40234 - 56874.73438 
56875.23828 - 56892.82422 
56890.65495 - 56909.82292 
56902.50000 - 56913.83203 



56980.07943 - 56995.10677 
57046.14063 - 57065.39063* 
NA 

57090.16797 - 57105.05469* 
57145.91406-57159.88672 
57250.90039 - 57264.17773* 

NA 

NA 

56903.37891 -56920.98828 
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3.3 Comments on individual candidates 

3.3.1 HD 23765 

This star is well isolated and the archival SuperWASP observations 
confirm the eclipse depth and show that only one very faint star re- 
mains as a potential contaminant. The host star is faint, for STEREO 
to observe such low amplitude variability, and detail on the shape of 
the eclipse is poorly defined but the SuperWASP data shows clear 
u-shaped eclipses, along with evidence of slight tidal distortion in 
the host. Given the spectral type of the host and the short period, 
there is the potential for this to be due to a very low mass compan- 
ion. The SuperWASP lightcurve is comparatively noisy and some 
cleaning was required in PERANSO, with slight indications that this 
might be partly due to very long period variability on a timescale 
greater than one year, although this would require more detailed ob- 
servations to confirm or refute. Undetrended STEREO observations 
do not show this but there are hints of sinusoidal variability with a 
period about seven times the orbital period, however STEREO/H1- 
1B data could not be included in this check for longer period vari- 
ability due to interference from several de-pointing events. 

3.3.2 HD 287039 

An apparently isolated star, there appear to be no candidates that 
could contaminate the lightcurve. The eclipses are fairly deep, be- 
ing 75 mmag, compared to others in the sample but also of short 
duration and the ingress and egress phases appear sharp. The T e ff 
of 5310 K estimated here results in a lower estimated radius for the 
host than all but one of the other stars in the sample. The eclipses 
are therefore consistent with a companion of 0.23 R@. There is no 
sign of tidal distortion in the host, despite the short period, and the 
companion may therefore be very low mass. 

3.3.3 HD 89849 

This star is well isolated with only a single very faint star near 
enough to potentially contaminate the lightcurve. The eclipses are 
box-like but with sufficient depth (54 mmag) and clear tidal distor- 
tion in the host that the companion is expected to be in the low mass 
stellar regime. The radius estimated for the companion is 0.23 R@, 
in agreement with this conclusion. However, this is dependent upon 
the accuracy of the estimated 1 Rq for the primary, which has no 
parallax but a T e ff of 5789 K. 

3.3.4 BD-07 3648 

For this star, the photometry used in calculating the effective tem- 
perature was taken from the NOMAD 1 catalogue (Zachari as et al.| 
|2004[ >. The star was initially suspected as being a low mass sys- 
tem due to the large parallax of 35.4 ± 19.6 mas jvan Altena, Lee| 
& Hoffleit 1995) combined with the faintness of the object. Un- 
like the other stars in this sample, however, the large eclipse depth 
indicates that the secondary has a similar size to the primary and 
therefore that the photometry might include a significant compo- 
nent of light from this star also. From the B — V colour we esti- 
mate an effective temperature of about 6000 K. Making use of the 
available 2MASS photometry we obtained three comparable val- 
ues of r e ff around 4000 K, therefore we finally adopted a Z^g of 
4074±9 1 8 K. From the transit depth, we know that the two stars 
have rather similar radii, meaning that the photometry might be af- 
fected by the secondary star, thereby introducing systematic errors 



in the adopted effective temperature. Taking into account the given 
parallax and the adopted T e g we estimate the stellar radius to be 
-~O.3/?0. However, contamination of the photometry by the light 
of the secondary star would result in an overestimation of the lu- 
minosity and temperature and, following from this, of the stellar 
radius, although the estimated radius does not seem unreasonable. 
Given the little information available and the complexity of the sys- 
tem, we are able just to give an indication of the spectral type of 
the two stars. Most likely the primary component is a late K-type 
star, while the secondary would then be an early M-type star. As a 
check, we attempted a fit of synthetic spectral energy distributions 
(SED), calculated with MARCS models jGustafsson et al.|2008) , 
to the available photometry. Using a single SED the best fit is ob- 
tained for a dwarf star with a r ef f of about 4250 K, which is then 
improved by adding a second SED for a dwarf star with a T e ff of 
about 3800 K, in good agreement with our conclusions. As a further 
check, we attempted to fit the photometry assuming a much hotter 
primary star with a T e g of 6000 K, as suggested by the B — V color. 
In this case, the fit of the infrared photometry was obtained only 
by assuming that the secondary star is a cool star with a radius of 
-~12^0, which is impossible given the observed transit depth, con- 
firming therefore our conclusions. This would be a rather challeng- 
ing system to confirm with spectroscopy or radial velocity studies. 
Nevertheless, as a nearby eclipsing binary of two very low mass 
stars it is more amenable to follow-up than the more distant sys- 
tems of this type. Secondary eclipses are just visible in the phase- 
folded lightcurve, showing no eccentricity. The secondary eclipse 
depth is hard to gauge, being visually of the order of 50 mmag deep 
whereas the scatter on the lightcurve is of the order of 0.35 mag. 

3.3.5 HD 75767 

This star is a known spectroscopic binary with an orbital period 
of 10.248042 ±0.000011 days, which has been determined by 
radial velocity measurements (Griffin 2006). The T e ff and abso- 
lute bolometric magnitude obtained from spectroscopic measure- 
ments in |Mishenina et al.| j2009| > are 5823 ±4 K and 4.60, re- 
spectively, which are in agreement with our photometric determi- 
nation of 5649 ±245 K and 4.55. The two stars are given to be 
spectral types GO and M3 and this agrees very well with the radii 
found from this study. The cause of a brightening in the phase- 
folded lightcurve around the primary eclipses is a systematic ef- 
fect due to the polynomial detrending, as illustrated by Figure [2] 
The radial velocity data clearly indicates that the eccentricity is 
very small at 0.0134 ± 0.0020 but secondary eclipses are not ob- 
served. HD 75767 is actually a quadruple system, with the spectro- 
scopic binary central pair orbited by a more distant M dwarf pair 
(Fuhrmann et al. 2005 1. This star has been known to be a spectro- 
scopic binary for a long time ( |Sanford|1931) but the eclipses have 
not been observed prior to STEREO/Hl-l, probably because of their 
shallow depth and short duration, the long period and also a lack of 
suitable reference stars for ground-based differential photometry. 
Figure[6]shows the radial velocity data collected by Griffin (1991) 
and |Fuhrmann et al.| p005l > for this star with the STEREO/Hl-l 
phase-folded lightcurve overlaid. There appears from Figure [6] to 
be a small discrepancy in the phase of the radial velocity mea- 
surements collected by Griffin ( 1991 1 and those of Fuhr mann et| 
[aT] ( |2005j > and it is tempting to speculate that this is due to the influ- 
ence of the circumbinary red dwarf binary reported by |Fuhrmann| 
|et al.| (2005). The earliest radial velocity measurements date back 
to the 1920's l|Sanford 1931) and it is not inconceivable that the 
circumbinary stars might perturb the system in this way on decadal 



8 K.T. Wraight, L. Fossati, Glenn J. White, A. J. Norton and D. Bewsher 



BD-D7 3648 P = 2.5265 DAYS 



HC 198044 P = 5.1520 DA V S 




I H.-.r-E 

HD 222891 P=1,5950 DAYS 



HD 213597 P = 2.4238 DA V S 



! H-' E 

HD 23765 P= 1.6865 DAYS STEREO 




I- l-IASE 

HD 287039 P=2-21 1 1 DAYS 



-v— v — 



HD 75767 P=10.2478 DAYS 



F i-ia-: E 

O aSB^-O 71 .1.0/81 "lAYS 




Figure 4. Phase-folded lightcurves of the low mass EB candidates found from the analysis of STEREO/HI-] data (see Section[3j and SuperWASP lightcurve 
(4th row, left column) of one of the candidates phase-folded on the same period as the STEREO/Hl-l counterpart (HD 23765, 3rd row, right column). The 
best-fitting lines shown bin 40 data points together with outliers more than 3 o from the mean of each bin culled. Two cycles are shown. 
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Figure 5. Close-up view of the eclipses of the low mass EB candidates. Top row: BD-07 3648, HD 198044, HD 205403, HD 213597 and HD 222891. Bottom 
row: HD 23765 (STEREO/HI- 1 ), HD 23765 (SuperWASP), HD 287039, HD 75767 and HD 89849. The best-fitting lines shown bin 40 data points together 
with outliers more than 3 a from the mean of each bin culled, except for HD 198044 and HD 205403 (both 20 points), HD 75767 (10 points), HD 222891 (30 
points) and HD 287039 (50 points). 



timescales. If this is the case, then transit timing variations would 
confirm it as well as providing orbits and masses for both binaries. 
It would also be worth monitoring for changes in eclipse depth in 
case the inclination of the system is similarly changing due to the 
combined motion of the two binaries. 



3.3.6 HD 198044 

This star has been reported by |Samus et al.| ( |2012[ > as a suspected 
variable with an amplitude of 50mmag, although the type of vari- 
ability is not specified. It appears to be very isolated, with no obvi- 
ous candidates to contaminate the lightcurve. Archival SuperWASP 
observations of this star have been made but are sparsely sampled in 
phase and the eclipses are not observed. This is the second-brightest 
star in the sample and also has the second-longest period, factors 
which make observing with SuperWASP difficult. The ingress and 
egress phases of the eclipses are very short indeed but there is some 
tidal distortion in the host, placing the companion in the low mass 
stellar regime. This star is one of the closest in the sample for those 
with a parallax available, at just 50 ± 1 parsecs. 



3.3.7 HD 205403 

This star is quite well isolated with only a single star almost two and 
a half arc-minutes away potentially contaminating the lightcurve, 
however this is almost exactly three magnitudes fainter in R. The 
eclipses appear to be total and the ingress and egress phases are 
very short but some tidal distortion is evident in the host, making 
the companion most likely to be low mass star. It is unclear whether 
or not secondary eclipses are seen for this star, as an extremely faint 
feature in the lightcurve could be easily confused with the signature 
caused by tidal distortion. 



3.3.8 HD 213597 

As detailed in Wraight et al. (2011), this well-isolated star shows 
eclipses due to a potential substellar mass companion. There ap- 
pears to be an indication of tidal distortion in the host but the very 
shallow, box-like transits suggest the companion to be very close 
to the brown dwarf regime. Interestingly, this is an early-type star 
with a slow projected rotational velocity, Vsin; = 40 ±4 km s 
(Gtebocki & Gnacihski 2003} . It is therefore a strong candidate 
to be a chemically peculiar star. There are, currently, only two 
known transiting planetary companions of early-type stars, WASP- 
33b ( |Smith et al.|[20TTT > and KOI-13 ( |Szab6 et al.|[20lT> , both of 
which are much fainter than this. A radius measurement has been 
made for this star, giving r = 2.039 ± 0.303/?q (Masana, Jordi & 
Ribas 2006). There is a hint of asymmetry in the transit, although 
this needs to be confirmed with more sensitive observations. Using 
the radius measurement above, we derive the radius of the compan- 
ion to be 0.35 ±O.31/?0. The host radius determined here, photo- 
metrically, is 1.96 ± O.O6i?0, from which we derive the radius of 
the companion to be O.36±O.12i? , which agrees with this. Given 
this apparent radius, a brown dwarf companion is unlikely, however 
this object is in close orbit of a star with T e ff of nearly 7000 K and 
could easily have become inflated as a result. Radial velocity mea- 
surements from Nordstroe m et al.| ( |1997| > show a reflex motion of 
just over 60 km s~ 1 , which suggests a low mass companion above 
the brown dwarf mass range. 



3.3.9 HD 222891 

A very isolated star, there is only one potential contaminant and this 
is almost four magnitudes fainter, so there is very good confidence 
in both the identity of the eclipsing variable and of the observed 
eclipse depth of 39 mmag. The eclipses are of short duration but, 
although the shape is not clear, both the matched filter model and a 
close-up of the eclipse (Figure[3] right plot, and Figure[5] top-right 
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Figure 6. Radial velocity data of HD 75767 collected from |Griffin| {199T) and |Fuhrmaim et al.|{2005) wi th the phase-folded lightcurve of STEREO/HI- 1 
overlaid in phase and set to km/s. The errors on the RV data points are the rms residuals given in Griffin 1 1991 1 and as specified for each data point from 
Fuhr mann et al.|{2005) . The larger errors are from photograhic observations originating in Sanford ( 1931). 



plot, respectively) suggest the eclipses are most likely total. More 
sensitive observations are required to confirm this and to search for 
secondary eclipses, which are not observed. There is a faint hint of 
tidal distortion in the host but, given its spectral type of F8 and the 
observed eclipse depth, the companion is expected to be very low 
mass. 



4 CONCLUSIONS 

The lack of bright eclipsing low mass, brown dwarf and exoplanet 
host stars is a major obstacle to detailed analysis of these interesting 
objects and to advancing the understanding of stellar formation and 
evolution. STEREO is an ideal platform for carrying out a survey 
for such objects and whilst the sample of nine presented here con- 
sists of companions likely to be too massive to be exoplanets, all are 
likely to be stars of very late spectral type (M or late K). In particu- 
lar, BD-07 3648 appears to be one of the brightest detached eclips- 
ing binaries where both stars are in the low mass regime. If follow- 
up observations show any of these candidates to have masses in the 
brown dwarf regime, however, it would be a big step forward for 
research into brown dwarfs and their atmospheres and other char- 
acteristics. Follow-up observations are required to ascertain the na- 
ture of all the candidates presented here, except for HD 75767 for 
which a wealth of data is available (Figure [§). There are relatively 
few detached low mass EB systems, with just six previously known 
double-lined detached eclipsing binaries where the companion has 
a radius less than 0.4 (Feiden & Chaboyer 2012). This sample 
of nine bright low mass EB candidates represents a considerable 
improvement in both the statistics and variety of low mass stars 
which can be summed up as follows: 

• With a brightness range of V = 6.59 mag to V = 11.3 mag, 
they are a sample suitable to be studied in great detail, refining 
parameters of both components with unprecedented precision. 

• A range of periods have been observed from 1.59495 days 
to 10.2478 days, that shows low mass stars can be found in a di- 
verse range of environments and provides a new benchmark against 



which to test models of stellar formation and magnetic field gener- 
ation. 

• The host stars for these EBs have a diverse range of tempera- 
tures, from 6925 K to 4074 K, providing an excellent comparison 
for models of stellar evolution. 

• Spectral types for the host stars range from F0 to late K, in- 
cluding an early-type star (HD 213597) and a star that itself appears 
to be low mass (BD-07 3648), reinforcing the observation that low 
mass stars can be found in almost any environment. 

• The lower limits of the radii of the companions (and also the 
primary for BD-07 3648) are determined to be in the range 0.22 Rq 
to 0.36 .Kq, providing 10 stars in bright eclipsing systems that are 
expected to be near to the boundary where stars become fully con- 
vective, challenging models of stellar structure to fully explain their 
properties. 
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